Scaling in Dynamic Susceptibility of Herbertsmithite and Heavy- Fermion Metals 
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We present a theory of the dynamic magnetic susceptibility of quantum spin liquid. The obtained 
results are in good agreement with experimental facts collected on herbertsmithite ZnCu3(OH)eCl2 
and on heavy-fermion metals, and allow us to predict a new scaling in magnetic fields in the dy- 
namic susceptibility. Under the application of strong magnetic fields quantum spin liquid becomes 
completely polarized. We show that this polarization can be viewed as a manifestation of gapped 
excitations when investigating the spin-lattice relaxation rate. 
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INTRODUCTION 

Landau Fermi liquid (LFL) theory is highly successful 
in the condensed matter physics. The key point of this 
theory is the existence of fermionic quasiparticlcs defin- 
ing the thermodynamic, relaxation and dynamic proper- 
ties of the material. However, strongly correlated Fermi 
systems encompassing a variety of systems that display 
behavior not easily understood within the Fermi liquid 
theory and called non-Fermi liquid (NFL) behavior. A 
paradigmatic example of the NFL behavior is represented 
by heavy-fermion (HF) metals, where a quantum phase 
transition (QPT) induces a transition between LFL and 
NFL QPT can be tuned by different parameters, 

such as the chemical composition, the pressure, and the 
magnetic field. Magnetic materials, in particular cop- 
per oxides and organic insulators, arc interesting sub- 
jects of study due to a quantum spin liquid (QSL) that 
can emerge when they approach QPT and are cooled to 
low temperature T. Exotic QSL is formed with such hy- 
pothetic particles as fermionic spinons carrying spin 1/2 
and no charge. A search for the materials is a challenge 
for condensed matter physics Q . The experimental stud- 
ies of herbertsmithite ZnCu3(OH)6Ci2 and the organic 
insulator EtMc 3 Sb[Pd(dmit)2]2 have discovered gaplcss 
excitations, analogous to excitations near the Fermi 
surface in HF metals, indicating that ZnCu3(OH)gCl2 
and EtMc3Sb[Pd(dmit)2]2 are the promising systems 
to investigate their QPTs and QSLs The ob- 

served behavior of the thermodynamic properties of 
ZnCu3(OH)eCl2 strongly resembles that in HF metals 
since a simple kagome lattice being strongly frustrated 
has a dispersionless topologically protected branch of the 



spectrum with zero excitation energy |14| - |l7| . This indi- 
cates that QSL formed by the ideal kagome lattice is 
located on the ordered side of the fermion condensation 
quantum phase transition (FCQPT) that is characterized 
by the presence of the spectrum with zero excitation en- 
ergy Q. This observation allows us to establish a close 
connection between QSL and HF metals whose HF sys- 
tems are located near FCQPT and, therefore, exhibiting 
an universal scaling behavior 0, [lj, G3] ■ As we are deal- 
ing with the real 3D compound ZnCu3(OH)gCl2 rather 
than with the ideal 2D kagome lattice, we have to bear 
in mind that the magnetic interactions and the presence 
of layers of nonmagnetic Zn 2+ ions separating magnetic 
kagome planes in the substance can shift the QSL from 
the initial point, positioning it in front of or behind FC- 
QPT. Therefore, the actual location has to be established 
by analyzing the experimental data. As a result, the lo- 
cation coincides with that of HF metals, and turns out 
to be at FCQPT [H, E3] , as it is shown in Fig. [Q Thus, 
FCQPT can be considered as QPT of ZnCu 3 (OH) 6 Cl 2 
QSL and both herbertsmithite and HF metals can be 
treated in the same framework, so that QSL is com- 
posed of fermions and these with zero charge and spin 
(7 = ±1/2 occupy the corresponding Fermi sphere with 
the Fermi momentum pp [2 , 14 , 17 1 . The ground state en- 
ergy E(n) is given by the Landau functional depending 
on the quasiparticle distribution function n CT (p), where 
p is the momentum. In spite of numerous experimental 
facts collected in measurements of inelastic neutron scat- 
tering spectrum and spin-lattice relaxation rates on her- 
bertsmithite, a theoretical understanding of how the dy- 
namical spin susceptibility of QSL behaves on approach- 
ing QPT and how it is affected by external parameters, 
such as the magnetic field, is still missing. 
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In this letter we employ the Landau transport equa- 
tion to construct the dynamical spin susceptibility. We 
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elucidate how the calculated susceptibility is affected by 
magnetic field and describe experimental facts collected 
on herbertsmithite and heavy-fermion metals. The ob- 
tained results are in good agreement with the facts and 
allow us to predict a new scaling emerging under the ap- 
plication of magnetic field in the dynamic susceptibility. 
Taking into account that QSL becomes completely polar- 
ized in strong magnetic fields, we show that this polar- 
ization can be seen as the presence of gapped excitations 
when investigating the spin-lattice relaxation rate. 



DYNAMIC SPIN SUSCEPTIBILITY OF 
QUANTUM SPIN LIQUID AND 
HEAVY-FERMION METALS 

To construct the dynamic spin susceptibility 
x(q,w,T) = x'(q,w,T) + i X "{q,uj,T) as a func- 
tion of momentum q, frequency ui and and temperature 
T, we use the model of homogeneous HF liquid located 
near FCQPT Q. To deal with the dynamic properties 
of Fermi systems, one can use the transport equation 
describing a slowly varying disturbance <5n CT (q, u>) of 
the quasiparticle distribution function no(p), and 
n = 6n + no. We consider the case when the disturbance 
is induced by the application of external magnetic field 
B = Bo + ABi(q, ui) with Bo being a static field and 
XBi a w-dependent field with A — > 0. As long as the 
transferred energy u < qpp/M* « u, where M* is 
the effective mass and \x is the chemical potential, the 
quasiparticle distribution function 7i(q, ui) satisfies the 
transport equation (l8j 

(qv p - Lu)Sn a - qv p ^ ^ /^(ppJ&V^pi) 



qvp^cr^ s (B + XBi) 



CT1P1 



(1) 



Here u B is the Bohr magneton and e p is the single- 
particle spectrum. We assume that B is finite but not so 
strong to lead to the full polarization of the correspond- 
ing quasiparticle band. In the field Bo, the two Fermi 
surfaces are displaced by opposite amounts, ±BqUb, and 
the magnetization M. = u_b(<5?t,+ — <5n_), where the two 
spin orientations with respect to the magnetic field are 
denoted by ±, and Sn± = ^2 p Sn±(p). The spin sus- 
ceptibility x is given by \ = 9M/dB\ B=Ba . In fact, the 
transport equation ([T]) is reduced to two equations which 
can be solved for each direction ± and allows one to cal- 
culate 8n± and the magnetization. The response to the 
application of A_Bi(q, uj) can be found by expanding the 
solution of Eq. (^Q) in a power series with respect to 
M*ui /qpF- As a result, we obtain the imaginary part of 
the spin susceptibility 



X"(q,w) =n% 



1 



2-rtq (1 



(2) 



where Ffi is the dimensionlcss spin antisymmetric quasi- 
particle interaction The interaction Fq is found to 
saturate at Fg ~ -0.8 [H, [H so that (1 + F a ) is pos- 
itive. It is seen from Eq. @ that the second term is 
an odd function of u>. Therefore, it does not contribute 
to the real part x' an d forms the imaginary part x" ■ 
Taking into account that at relatively high frequencies 
lu > qpp/M* <§C u in the hydrodynamic approximation 
X 1 oc \/uj 2 [2lj . we conclude that the equation 



x(q,w) 



I'B 



M*p F 



7T 2 (1 + F a ) 1 



M*LJ 



(3) 



produces the simple approximation for the susceptibility 
X and satisfies the Kramcrs-Kronig relation connecting 
the real and imaginary parts of x- 

To understand how can x" arL d x given by Eqs. ([2]) and 
((3]), respectively, depend on temperature T and magnetic 
field B, we recall that near FCQPT point the effective 
mass M* depends on T and B, and is given by the Lan- 
dau equation (LE) [H, (UJ . The interaction function F of 
LE is completely defined by the fact that the system has 
to be at FCQPT. The sole role of F is to bring the system 
to FCQPT, where the Fermi surface alters its topo logy 
so that M* acquires T and B dependencies 0, Uj, l24j . 
At FCQPT, LE can be solved analytically: At B = 0, 
the effective mass depends on T 



M*(T) ~ a T T 



-2/3 



(4) 



At finite T, the application of magnetic field B drives the 
system to the LFL region with 



M*(B) ~ a B B- 2/3 . 



(5) 



Here ax and as are constants. At finite B and T near 
FCQPT, the solutions of LE can be well approximated 
by a simple universal interpolating function. The inter- 
polation occurs between the LFL (M*(T) oc const) and 
NFL (M*(T) oc T~ 2/3 ) regions. It is convenient to intro- 
duce the normalized effective mass and the normal- 
ized temperature Tn dividing the effective mass M* by 
its maximal values, M^ ax , and temperature T by T max 
at which the maximum occurs. The normalized effective 
mass = M* /M^ ax as a function of the normalized 
temperature y = Tn = T/T max is given by the interpo- 
lating function Q 



M* N (y) » c 



l + c lV 2 



1 



c 2 y 



8/3 ' 



(6) 



Here Co = (1 + C2)/(l + c%), c\ and C2 are fitting pa- 
rameters, making M^(y = 1) = 1. Magnetic field B 
enters LE only in the combination fisB/ksT, making 
fcsTmax — (J-bB where ks is the Boltzmann constant 
0, Ull- Thus, in the presence of magnetic fields the vari- 
able y becomes y = T/T max ~ fc^T '/ ' [IbB. Since the vari- 
ables T and B enter symmetrically Eq. ^ is valid for 
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FIG. 1: (color online). T — B phase diagram of QSL and 
HF liquid. The vertical and horizontal arrows, crossing the 
transition region depicted by the thick lines, show LFL-NFL 
and NFL-LFL transitions at fixed B and T, respectively. At 
temperatures T < T* and magnetic field B < B* shown by 
the dash-dot arrows the effective mass M* ~ const and the 
system in the LFL region. The dash line continuing the thick 
line represents the transition region provided the system were 
located at FCQPT shown by the arrow. 



y = I-IbB /ksT '. Now we construct the schematic T — B 
phase diagram of QSL and HF liquid reported in Fig. [T] 
At T = and B = the system can exactly be located 
at the FCQPT point without tuning with both T* and 
B* arc zero. It can also be shifted from the FCQPT 
point by doping, pressure etc. In that case T* and B* 
become finite so that at T < T* and B < B* the effec- 
tive mass M* ~ const. As seen from Fig. Q] at T ~ T* 
and B ~ B* the transition region exhibits a kink, since 
M* is no longer constant at rising B and T. Magnetic 
field B and temperature T play the role of the control 
parameters, driving it from the NFL to LFL regions as 
shown by the vertical and horizontal arrows. At fixed 
B and increasing T the system transits along the verti- 
cal arrow from the LFL region to NFL one crossing the 
transition region. On the contrary, at fixed T increasing 
B drives the system along the horizontal arrow from the 
NFL region to LFL one. 
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SCALING BEHAVIOR OF THE DYNAMIC 
SUSCEPTIBILITY 

To elucidate a scaling behavior of x, we employ Eq. (0| 
to describe the temperature dependence of x- It follows 
from Eqs. © and @ that 



T 2/3 x(T,lu) 



a i 



1 + ia 2 E 



(7) 



Here a\ and a 2 arc constants absorbing irrelevant values 
and E = w/(fcsT) 2 / 3 . As a result, the imaginary part 



FIG. 2: (color online). The function (T 2/3 x")n plotted 
against the unitless ratio En = oj/((fcsT) 2//3 i5 max ). The data 
extracted from measurements on ZnCu3(OH)eCl2 obtained 
for 0.077 < T < 42 K 0, Panel A, and on the HF metal 
Ceo.925La .o75Ru2Si 2 obtained for 2.5 < T < 80 K at Qi [25l |. 
Panel B and C, collapse onto a single curve. The solid curves 
are fits with the function given by Eq. ©. 
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x"(T,u>) satisfies the equation 
T 2 / 3 X"(7» 



a 3 E 



1 + a 4 E 2 ' 



(8) 



where a 3 and a 4 are constants. It is seen from Eq. ([8]) 
that T 2 / 3 X "(T,uj) has a maximum (T 2 / 3 x"(T, w)) max at 
some .Emax and depends on the only variable E. Equa- 
tion (EJ) confirms the scaling behavior of ^"T 6e experi- 
mentally established in Ref. Q. As it was done for the 
effective mass when constructing ([6]) , we introduce the di- 
mensionless function (T 2 ^ X ")n = T 2 / 3 x7(T 2 / 3 x")max 
and the dimensionless variable En = E/E max , and Eq. 
(El) is modified to read 



{T 2/ \") 



biE. 



N 



N 



l + b 9 E 



2 ' 

N 



(9) 



with bi and b 2 are fitting parameters which are to ad- 
just the function on the right-hand side of Eq. © to 
reach its maximum value 1 at E n = 1. We predict that 
if measurements of X " are taken at fixed T as a func- 
tion of B, then taking into account Eq. ([5]), we again 
obtain that the function B 2 / 3 x"(E) exhibits the scaling 
behavior with E = w/(/ig.B 2 / 3 ). If the system is placed 
at FCQPT, the scaling described above is valid down to 
lowest temperatures. When the system is shifted from 
FCQPT then T* and B* arc finite and the scaling is vi- 
olated in the LFL region and recovered in the NFL one 
at T > T* or B > B* as seen from Fig. [Q 

In Fig. [2] consistent with Eq. , the scaling of the 
normalized dynamic susceptibility (T 2 / 3 x")n extracted 
from the inelastic neutron scattering spectrum of both 
herbertsmithite ML Panel A, and Ceo.925Lao.o75Ru2Si2, 
Panel B and C, [251] is displayed. The scaled data col- 
lapse fairly well onto a single curve over almost three 
decades of En ■ It is seen that our calculations shown by 
the solid curves are overall in good agreement with the 
experimental facts. We note that, as seen from Fig. [21 
Panel C, the data taken at 2.5 K are at variance with 
the scaling behavior shown by the solid curve. We sug- 
gest that Ceo.925Lao.075Ru.2Si2 is slightly shifted from 
FCQPT as shown in Fig. [1] and at sufficiently low tem- 
peratures T < T* the scaling is violated [ijj [26| , while 
ZnCu3(OH)6Cl2 is near that point. Some remarks on a 
role of both the disorder and the anisotropy are in or- 
der. The anisotropy is supposed to be related to the 
Dzyaloshinskii-Moriya interaction, exchange anisotropy, 
or out-of-plane impurities. Measurements of the suscepti- 
bility on the single crystal of herbertsmithite have shown 
that it closely follows that measured on a powder sample 
Q. At low temperatures T < 70 K, the single- crystal 
data do not show magnetic anisotropy These con- 
firm that the stoichiometry, disorder and anisotropy do 
not contribute significantly to the results at relatively 
low temperatures. Moreover, the scaling behavior of the 
thermodynamic functions of herbertsmithite is the in- 
trinsic feature and has nothing to do with the impurities 



_17j . These observations are in agreement with a general 
consideration of scaling behavior of HF metals Q . 



SPIN-LATTICE RELAXATION RATE OF 
QUANTUM SPIN LIQUID 

Consider the effect of B on the spin-lattice relaxation 
rate \jT{T determined by X " given by Eq. ((5J) 



X"(q,< 



-| u _>ooc(M*) 2 , (10) 



where Aq is the hyperfine coupling constant of the muon 
(or nuclei) with the spin excitations at wave vector q 
[H, [27], [28| . Figure [3] and the inset display the normalized 
(1/T\T)n and the normalized longitudinal magnetore- 
sistance pn at fixed temperature versus the normalized 
magnetic field Bn- It is seen from Fig. [3] that the mag- 
netic field progressively reduces \jT{F and the longitu- 
dinal magnetoresistance (LMR) , and these as a function 
of B possess an inflection point at B — Bi n f shown b 
the arrow. The normalized LMR obeys the equation 

„ /o % _ p( b n) - A) _ ( 1 \ _/ M *yi 

Pn{B N ) - - 77T7= - (M N ) , 

Pinf J n 



2] 



(11) 



where po is the residual resistance, pi n f is LMR taken 
at the inflection point, p is LMR, and Bn = B/Bi n f. 
We normalize (1/TiT) and LMR by their values at the 
inflection point, and the magnetic field is normalized 
by B in f. In accordance with the phase diagram [TJ at 
B > B in f, as seen from Fig. [31 QSL enters the LFL 
region with independence of the effective mass defined 
by Eq. ©. It follows from Eqs. ([TDJ and (fTTJ) that 
(1/TiT)jv = Pn = (M N ) 2 where {M N ) 2 is defined by Eq. 
([6]) which shows that different strongly correlated Fermi 
systems are to exhibit the same scaling of {M N ) 2 . It is 
seen from Fig. [3] and from the inset, that YbCu5_ x Au x , 
herbertsmithite ZnCu3(OH)6Cl2 and YbRIi2Si2 demon- 
strate the similar behavior of (M N ) 2 resulting in the scal- 
ing of LMR and l/T x T. Thus, Eqs. ©, ^ and {TT) de- 
termine the close relationship existing between the quite 
different dynamic properties and different strongly corre- 
lated Fermi systems such as QSL and HF metals, reveal- 
ing their scaling behavior at FCQPT. 

We note that one may be confused when applying Eq. 
(fTU| to describe (1/TiT) in strong magnetic fields. In 
that case both QSL and HF metals become fully polar- 
ized due to Zeeman splitting 14, 27, 28, |32|. As a re- 
sult, one subband becomes empty, while the energy ef 
of spinons at the Fermi surface of the other subband lies 
below the chemical potential p formed by the magnetic 
field B . It follows from Eq. d) that X " = and Eq. 
(fTT?ll is not valid. The difference 6 = p — Ef can be 
viewed as a gap that makes 1/TiT oc cxp — (5/ksT). At 
temperatures fc^T ~ S, the subbands arc populated by 



5 




FIG. 3: (color online). Magnetic field dependence of normal- 
ized (see text for details) muon spin-lattice relaxation rate 
(1 /TiT) n extracted from measurements on YbCu4.4Auo.6 [29l ] 
and ZnCu3(OH)6Ci2 [13] along with the normalized longitu- 
dinal magnetoresistance pjv versus normalized magnetic field 
Bn- Our calculations are shown by the solid line. The ar- 
rows indicate the inflection points. Inset: pM versus Bn, Pn is 
extracted from measurements on YbRli2Si2 at different tem- 
peratures [U listed in the legend. The solid curve represents 
our calculations. 

spinons and the validity of Eq. (JTUJ) is restored. Thus, 
5 can be interpreted as the presence of gapped excita- 
tions. On the other hand, if there were the gapped ex- 
citations, then the heat capacity demonstrates the ex- 
ponential decay rather than a linear T-dcpcndcnce at 
low temperatures. Analysis based on experimental data 
shows the presence of linear T-dcpcndence even under 
the application of high magnetic fields [3] , while recent 
measurements on ZnCu3(OH)gCl2 of \jT\T suggest the 
gapped excitations (33|. To clarify whether the gapped 
excitations would occur in ZnCu3(OH)gCl2, an accurate 
experimental measurement in magnetic fields of the low 
temperature heat capacity is necessary. 



CONCLUSIONS 

We have presented a theory of the dynamic magnetic 
susceptibility of quantum spin liquid, and elucidated 
how the calculated susceptibility is affected by magnetic 
field and describe experimental facts collected on her- 
bertsmithite and heavy-fermion metals. The obtained 
results are in good agreement with experimental facts 
collected on both herbertsmithite ZnCu3(OH)gCl2 and 
on heavy-fermion metals, and allow us to conclude that 
the dynamic magnetic susceptibility of herbertsmithite 
is similar to that of heavy-fermion metals. Thus, her- 
bertsmithite can be viewed as a new type of strongly 



correlated electrical insulator that possesses properties 
of heavy-fermion metals with one exception: it resists 
the flow of electric charge. We have also predicted a 
new scaling in magnetic fields in the dynamic suscepti- 
bility emerging under the application of magnetic field. 
Taking into account that under the application of strong 
magnetic fields quantum spin liquid becomes completely 
polarized, we have shown that this polarization can be 
viewed as a manifestation of gapped excitations when in- 
vestigating the spin-lattice relaxation rate. 
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